Theta-Defensins to Counter COVID-19 as Furin Inhibitors: In Silico Efficiency Prediction and Novel Compound Design

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was characterized as a pandemic by the World Health Organization (WHO) in Dec. 2019. SARS-CoV-2 binds to the cell membrane through spike proteins on its surface and infects the cell. Furin, a host-cell enzyme, possesses a binding site for the spike protein. Thus, molecules that block furin could potentially be a therapeutic solution. Defensins are antimicrobial peptides that can hypothetically inhibit furin because of their arginine-rich structure. Theta-defensins, a subclass of defensins, have attracted attention as drug candidates due to their small size, unique structure, and involvement in several defense mechanisms. Theta-defensins could be a potential treatment for COVID-19 through furin inhibition and an anti-inflammatory mechanism. Note that inflammatory events are a significant and deadly condition that could happen at the later stages of COVID-19 infection. Here, the potential of theta-defensins against SARS-CoV-2 infection was investigated through in silico approaches. Based on docking analysis results, theta-defensins can function as furin inhibitors. Additionally, a novel candidate peptide against COVID-19 with optimal properties regarding antigenicity, stability, electrostatic potential, and binding strength was proposed. Further in vitro/in vivo investigations could verify the efficiency of the designed novel peptide.


Introduction
The new coronavirus, called "severe acute respiratory syndrome coronavirus 2" (SARS-CoV-2), was detected in China in Dec. 2019 for the first time as the agent of the coronavirus disease 2019 (COVID-19) outbreak [1], which was announced as a pandemic by the World Health Organization (WHO) in March 2020. It has caused more than 276 million confirmed cases and more than five million deaths globally as of December 2021 [2]. Numerous clinical trials were started to investigate the efficacy and safety of various drugs such as hydroxychloroquine, lopinavir/ritonavir, corticosteroids, and remdesivir for this disease, and guidelines are being updated continually [3]. Combinations of drugs with different mechanisms might be more effective [4]. Fortunately, several vaccines have been developed and are being used, but antiviral drugs are still needed to save the lives of infected people. Moreover, various mutations on the SARS-CoV-2 genome indicate the need for drugs besides vaccines. Treatment approaches are divided into two broad categories according to the targets: (1) inhibition of enzymes involved in viral pathogenesis and (2) inhibition of lung injury that happens due to inflammation by immune system modulation [5].
Coronaviruses have spike glycoproteins on their surface, which bind to the cell membrane and infect the cell [6]. A large ectodomain, a single-pass transmembrane anchor, and a short C-terminal intracellular tail are the components found on the spike glycoprotein [7]. The ectodomain consists of a receptor-binding unit, called S1, and a membrane-fusion unit, named S2 [8]. For the virus entry into host cells, binding of S1 to a specific cell surface receptor through its receptor-binding domain (RBD) is the first step, followed by S2 fusion with the host cell, which finally leads to infusion of the viral genome into host cells [9].
Specific cell enzymes such as furin, which has a site on the spike protein, can facilitate the entry process of the viral genome [8]. Furin is a host-cell enzyme classified as a family member of proprotein convertases. A total of nine members of this family have been identified in humans [10]. Furin is a proprotein convertase subtilisin/Kexin type 3, also termed PACE (paired basic amino acid cleaving enzyme), because it cleaves basic amino acid motifs. Furin regulates the activity of various mammalian, bacterial, and viral proteins. Therefore, it can be considered a target for treating different infectious and noninfectious diseases [11]. This host-cell enzyme is regarded as necessary in the viral maturation process and seems to be involved in SARS-CoV-2 pathogenesis and probably viral transfer among humans [12]. Thus, furin protein blockade is suggested as a therapeutic strategy against SARS-COV-2 [13].
Several peptidyl and nonpeptidyl furin inhibitors have been developed in the last decade. Some have shown promising results in mouse models [14], but only a few have entered human clinical trials for cancer treatment [11,15]. The promising results of peptide-based furin inhibitors in treating viral infections have been reviewed [16].
Defensins are arginine-rich antimicrobial peptides, which can hypothetically inhibit furin since it was proven that poly-arginine-derived molecules could inhibit the activity of furin protein [17]. Plants, animals, and fungi express defensins as a defense tool [18,19], regarded as essential for the host response to infections. Defensins have antimicrobial effects against Gram-positive and Gram-negative bacteria, fungi, parasitic protozoa, and enveloped and nonenveloped viruses [18,20]. They target several different stages of the virus life cycle [21], such as viral replication [22], and can bind to viral or host proteins [23,24]. Defensin can also modulate innate and adaptive immune responses [25].
Defensins are usually cationic with multiple cysteines in their structure [18], which secure a beta-sheet core structure through three conserved intramolecular disulfide bonds [26]. There are glycine and positively charged amino acids (arginine and lysine) in the defensin structure, as well as aliphatic hydrophobic residues that form a hydrophobic core [27].
Three families of mammalian defensins are identified based on their structure, including alpha-, beta-, and thetadefensins (TDs) [28]. TDs contain several arginine residues and cyclic peptides in their structure. They are expressed in some animal species, such as Old World monkeys, and cannot be found in human cells because of a premature stop codon [29]. Synthetic humanized TDs could have a therapeutic effect on bacterial and viral lung infections [30] and have anti-infective and anti-inflammatory properties [31]. These peptides act against viruses (such as HIV and influenza) by various mechanisms such as direct interaction with the virus [32,33], protease inhibition [34], and suppression of proinflammatory responses [35]. A rhesus type of TD, Rhesus theta-defensin-1 (RTD-1), protected mice against SARS infection, a possible mechanism of which was the immunomodulatory activity of RTD-1 [36]. TDs have attracted attention as anti-infective drug candidates due to their small size, unique structures, and involvement in several defense mechanisms [29].
Bioinformatics tools help study different compounds and predict their features in less time and lower costs than the lab experiments, giving the insight to find or design more efficient structures [37,38]. Computational approaches could be of great value given the importance and urgency of finding therapeutics molecules for COVID-19. This study investigates the potential of TD against SARS-CoV-2 through an in silico analysis of the TD structure and delves into the nature of TD-furin interaction by docking approaches. Finally, a candidate novel peptide against COVID-19 with furin inhibition and anti-inflammatory activity would be proposed.
To find any peptide-reactive furin interactions, querying the amino acid sequence of furin was also conducted in PepBDB. PepBDB is a database containing information about biological peptide-mediated protein interactions of the peptides in the PDB with a length of up to 50 residues.
2.2. Sequence Comparison. Structural alignment of TD sequences was performed by the T-Coffee expresso algorithm at http://tcoffee.org/ [43] for sequence comparison. The duplicated sequences were removed before alignment by the duplicate finder java standalone application. The alignment file was visualized and explored by the Alvis alignment visualizer [44].

TD Consensus
Patterns. The conserved patterns on TD sequences were discovered using the PRATT algorithm [45] as provided at https://web.expasy.org/pratt/. The PRATT tool enables finding conserved patterns in a group of protein sequences.
2.4. Antigenicity. The antigenicity of each defensin was assessed by the VaxiJen server [46] at http://www.ddgpharmfac.net/. VaxiJen helps to search for protective antigens through an alignment-independent method, which performs the prediction about the physicochemical properties of proteins.
2.5. Electrostatic Potential. The electrostatic potential of each defensin peptide and furin was measured based on the electrostatic surface grid calculated according to Coulomb's law [47]. First, the molecular surfaces were displayed by the 2 Computational and Mathematical Methods in Medicine UCSF Chimera program [48], which could be used for molecular modeling, visualization, and analysis. The surface-enclosed residues were used to calculate the electrostatic potential. The distance-dependent dielectric was selected to vary in proportion to the distance from each charge. The dielectric constant was set at 4. The distance from the surface vertex was 1.4; the protonation state was based on a local H-bonding environment. The surface coloring was performed by the "Coulombic Surface Coloring" function of the UCSF Chimera program [48]. The molecular surface coloring is performed in this method by the potential values. It can be used for structures with or without explicit hydrogens and could produce a grid of potential values.
2.6. Investigating the Furin Properties. The properties of furin active sites were assessed by the CAVER Analyst version 2 standalone software [49]. The approach was initiated by defining the active sites extracted from the databases. The volume data was based on the Voss Volume Voxelator (3V) algorithm [50], as provided at http://3vee.molmovdb.org/. The 3V encompasses several programs for estimating volumes in protein files.
The properties of the furin cleft were assessed by the CavityPlus server [51] at http://www.pkumdl.cn/. It helps to detect and analyze the functions of protein cavities using 3D structural information of proteins.
2.7. Docking. The binding properties of furin and the available TDs were defined by applying a hybrid approach strategy (template-based modeling into free docking with default parameters) by HDOCK [52] (http://hdock.phys .hust.edu.cn/).
A peptide-based docking approach was employed on all peptide sequences separately by GalaxyPepDock, available at http://galaxy.seoklab.org. The server employs a combination of information on similar interactions found in the structural databases and model building through energy-based optimization [53].
2.8. Network View of Peptide-Furin Interactions. The networks of the interaction of each available TD with furin were extracted from the docked structures. First, the TD peptide structure was selected within the protein complex; then, the network of interactions between the selected residues and other residues in the complex was extracted by using RINanalyzer [54], as embedded in Cytoscape 3.7.2 [55]. The RINanalyzer plugin of Cytoscape allows creating an interaction network of amino acid residues directly from the structure by connecting the UCSF Chimera with Cytoscape. The residue interaction networks included hydrogen bonds, combined edges, contacts, clashes, backbone connectivity, and distances. The extracted networks were visualized and explored by using the Cytoscape software tool.  Computational and Mathematical Methods in Medicine protein structure) [56]. The approach was made by the RINspector application [57] of Cytoscape. A z-score was calculated for each node based on modifying the average shortest path length upon removal of each node. The central residues, which are the residues that play an essential role in the communication within the interaction network, were assigned by a significant z-score (z-score > 2).
2.10. Nature-Inspired Peptide Design 2.10.1. Construct Properties and Modeling. The tertiary structure of the construct was built by the de novo method at https://bioserv.rpbs.univ-paris-diderot.fr/services/PEP-FOLD/ by the PEP-FOLD server [58] and the Modeller software [59] using the graphical user interface available in Chimera. The structure prediction was performed using a homology modeling approach. The quality of the predicted structures was evaluated by MolProbity [60] available at http://molprobity.biochem.duke.edu/.
The Disulfide by Design, provided at http://cptweb.cpt .wayne.edu/, was employed to examine the disulfide bonds in the predicted structure.
The construct dynamics were determined and visualized by sampling conformations using the normal analysis mode of the Dynamut server [61] (http://biosig.unimelb.edu.au/ dynamut/). The C-alpha force field was selected for standard mode analysis calculations, which was derived from fitting to the Amber94 all-atom potential.
The designed construct was subjected to the aforementioned evaluations, including antigenicity, stability, electrostatic potential, and binding strength.

Data Mining.
Due to the presence of several solved structures in PDB for the 794-amino acid reference sequence of furin (UniProtKB ID: P09955), the amino acid sequence of furin was obtained using the archived cluster of structures (PDBFlex database). The results included a range of clusters at a similarity level of 49.01 (cluster master: 6F9M) to 99.36 (cluster master: 6HZA). 6HZA was comprised of 43 members, among which the best resolution was related to the X-ray diffraction solution of furin under the accession number of 5JXG and resolution of 1.80 Å.

Furin
Properties. The reactive residues, active site, and active cleft of furin were obtained using databases and literature. The available UniProtKB data on the furin reference sequence is summarized in Table S1 of Supplementary File 1. D153, H194, and S368 were introduced as active sites based on the charge relay system (Prosite rule). Furthermore, the interactive residues in the furin-peptide complexes,   Table S2). These data provided a reference to confirm the docking results and helped define the interactive residues in the docking process.
3.2.1. Furin Cleft. All cavities on the structure surface were identified to measure the volume and area of the furin cleft. At least nine cavities were recognized on the surface of furin (PDB entry 5JXG, Figure S1). A large spherical cavity with a volume of 1968 Å 3 , an area of 1371 Å 2 , and sphericity of 0.8 presents the active cleft of the enzyme. The route to the active site also passed through this relatively large cavity ( Figure 1). The maximal pKd of the binding site was 10.13, which revealed the druggability of this cavity. The electrostatic measurement of the surface area showed the extremely negative electrostatic potential of the active cleft (-28) ( Figure 2). As shown in Figure 2, the active cavity and its surroundings comprised a negative electrostatic potential, whereas the other areas showed a medium potential. The high potentials were scattered through small surfaces (green regions in Figure 2).

TDs
3.3.1. Sequence Analysis. All defensins used in the present study were short (approximately 18-residue long) peptides. The existence of several arginines and cysteine residues in the sequences of all defensins is mentionable. The arginine and cysteines are both the fundamental elements that appear in the sequence pattern. All sequences obey a consensus pattern, namely, R-C- x is any amino acid; this pattern served as a template for designing a novel construct.
The stability and antigenicity of each peptide were based on the physicochemical properties and amino acid composition of the sequences. Table 1 shows the sequence stability and antigenicity of each peptide in which instability indices lower than 40 are considered stable, and a lower number offers more stability. Thus, the most stable TD was 2ATG followed by 2LZ1, while other peptides were indicated to be unstable. However, these two peptides were presented with high antigenicity. The antigenicity threshold is 0.4, and the fewer antigenic peptides were 2M2H, 2M2X, 2M2G, 1HVZ, and 2LYE, respectively. Our designed construct had a low antigenicity with close-to-threshold stability.

Structural Exploration.
The list of TDs used in this study is as follows. The PDB IDs of TDs are 1HVZ, 2LYE, 2LZI, 2M1P, 2M2G, 2M2H, 2M2S, 2M2X, 2M2Y, 2M77, 2M78, and 2M79. These theta-defensins include both natural and synthetic molecules. Since all these structures are determined by nuclear magnetic resonance spectroscopy and contain several trajectories, the model for further analyses was selected based on the full validation report deposited in PDB. The medoid structure was selected. The   [62,63], is the model that is the most similar to all other models.

The Convexity Indices.
A convexity index was assigned for each residue to define the protruded side chains of the amino acids of TDs. The most protruded side chains were attributed to arginine (Supplementary Data 1, Table S3; Supplementary Figure S2).

The Electrostatic Potential.
A positive electrostatic complement peptide would engage the highly negative furin cleft. To account for this, along with calculating the furin electrostatic potential, similar values were calculated for each peptide. The results showed that all defensins were highly favorable in electrostatic potential (Supplementary Figure S3). This would explain the binding property of these peptides to the furin cleft. Moreover, the peptides' calculated volume and surface area revealed the appropriate size of these peptides to occupy the furin cleft (Table 2).
3.3.5. Furin-Peptide Docking. Two docking approaches were performed to delve into the nature of theta-defensin-furin interaction. The first was based on the structures, and the second was on the peptide sequence and structure of furin. The structure-based method resulted in ten different orientations for each peptide structure and furin structure.
Considering the docking scores, it can be observed that the minimal energies were not necessarily related to the correct orientation (Supplementary Data 1, Table S4). Our primary criterion defines the correct orientation (the engagement of the peptide structure in the furin active cleft). The minimal energetic interaction of correct orientation belonged to the 2M1P structure (-271.05).

Essential Residues in Protein Complexes.
A centrality analysis outlined the essential residues involved in the peptide-furin complex. The z-score values (significant zscore: >2) of centrality reveal a notable role for arginine residues in each complex (Table 3 and Supplementary Figure S4). Instead, the cysteines were not attributed to significant z-scores. This emphasized the critical role of arginine residues in the peptide-protein interaction. In contrast, cysteines are critical in maintaining the structure of the peptides.

Rational
Designing of the Construct. The data obtained from the previous sections provide a guideline for designing a novel peptide, which will be referred to as "construct" hereafter. The criteria for designing the construct are achieving a peptide with higher stability, lower antigenicity, and higher electrostatic potential than natural TDs. The sequence pattern extracted from the library of TDs provided a scaffold for rational designing of the novel construct. In the designing plan, the construct sequence was   The sequence of the designed construct was also included in the alignment. Each line represents a peptide sequence that passed through the amino acid matrix introducing the peptide residue of the related site. Each sequence in the alignment block can be identified by following each line. When the bundles stack on an amino acid, a conserved site is introduced. The arrangement of amino acids was based on physicochemical properties specified at the left. At the bottom, the sequence logo also shows the conservancy of each side; the size of fonts is proportional to the entropy in the alignment column. Below the sequence logo, the sequence of the designed construct is presented by color codes.

Computational and Mathematical Methods in Medicine
The pattern remained constant, and the variable residues were selected from the options, which were limited to the existing residues identified in the sequence alignment ( Figure 3).
The proper residues were selected based on centrality z -scores; the residues of the significant z-score were selected within other options. Moreover, the scarcity of the residues in the antibody-antigen interface was followed to lower the antigenicity of the construct.
3.6. Construct Structure. The de novo modeling method and homology modeling approach revealed slightly different structures, yet the homology modeling approach yielded a more satisfactory structure (Supplementary Data 1, Table S5). In the homology modeling approach, two structures of 2LYE and 5NIZ served as templates with 88.9 percent sequence identity. The resulting modeled structures comprised the ZDOPE score of 0.62 to 1.7. The best structure (with the lower ZDOPE value) was selected for further investigations. In the modeled structure, the three expected disulfide bonds are appropriately located with χ 3 angles of 101.32, 116.66, and 117.09; the energy of 0.73, 2.20, and 1.48 kcal/mol; and Σ B-factor of 143.08, 158.14, and 429.42 (the χ 3 angle is typically -87°or +97°± 30). The predicted structure of the construct with its residue side chains and disulfide bonds is presented in Figure 4.
In Figure 2, the right panel also indicates the too high electrostatic potential of the construct (+10).
The convexity index list of atoms in the construct structure showed the high protrusion value of arginine side chains ( Figure 5 and Supplementary Data 1, Table S3); the highest sum of all-atom convexity index belonged to Arg10 (41.95), and the most protruded atoms were hydrogens (6.19).
To explore the properties of construct-furin interaction, the peptide construct was docked to furin. The minimal energy within the top ten predictions represented the correct orientation of the construct ( Table 4). The template for docking was furin in the complex with a competitive inhibitor (meta-guanidinomethyl-phenyl acetyl-Arg-Val-Arg-(amidomethyl) benzamidine) under the PDB accession number of 4OMC [64].
The interaction network of the construct-furin complex showed that the selected residues had significant z-scores ( Figure 6, Table 3).
To examine the flexibility of the construct, the molecular fluctuation was predicted. The vectors, trajectories, and fluctuations were merged to represent more flexible regions and the movement vectors (Figure 7). The very ends of the molecule comprised the highest level of fluctuations; the trajectories in Figure 7 evidenced the highly fluctuating (therefore, more flexible) arginine residues. Additionally, arginine and

Discussion
Given the role of furin in SARS-CoV-2 pathogenicity and the fact that inhibition of host enzymes could be a better approach than viral enzymes, furin is believed to be an important therapeutic target in COVID-19 management. Some findings support this theory. For instance, genetic variations of furin cause different binding affinities toward viral spike glycoprotein and result in various infectivity rates in humans [65]. Additionally, the plasma level of furin was found to increase in patients with coronary artery disease and COVID-19, which may be the reason for their poor clinical prognosis [66]. TDs are considered furin inhibitors and thus could be employed in the treatment of some viral infections, namely, COVID-19. In this study, the structure of TDs and their interaction with furin were examined in detail to verify the abovementioned notion and help design a similar construct with more beneficial features.
The low electrostatic potential of the furin active cleft and high electrostatic potential of the TDs, as well as the appropriate size of these peptides for insertion into the furin cleft, are evidence of a strong interaction between furin and TDs, which could lead to furin inhibition by TDs.
The observed pattern of TDs provided a clue for modifying their structure to obtain an improved TD with better binding property, which can engage with the furin cleft. Mining the data of the available structures provided clues for designing a novel peptide that would better occupy the furin cleft and interact more strongly besides possessing less immunogenicity, more stability, and more flexibility. The original template for designing the construct was the basic sequence pattern consensus within the available TDs. The results showed that peptide modifications led to a novel construct with more stability and hydrophilicity as well as a smaller size than natural TD. The physicochemical properties of the novel TD suggested it as an appropriate furin inhibitor or, in other words, a more robust binder.
Arginine residues have an essential role in TD-furin interaction. This result corresponds to the hypothesis that defensins, as arginine-rich molecules, can inhibit furin. On the other hand, cysteines are critical in maintaining the structure of TD. The arginine residues are protruded from two side loops that connect a central beta-ladder. The three cysteine pairs are coupled face to face in the beta-ladder to form three disulfide bonds. Arginine and isoleucine residue fluctuations cause the flexibility of the molecule. The nonpolar amino acids help the hydrophilic-lipophilic balance.
According to Table 1, most of the TD peptides represented antigenicity. Some structures, such as the 2M2H structure, had low antigenicity but were not stable enough. The minimal energy of interaction of correct orientation is related to the 2M1P structure (-271.05). Nevertheless, its antigenic properties and stability were not suitable. Therefore, an optimum state should be considered. The designed peptide (construct) had higher stability, lower antigenicity, and higher electrostatic potential. Generally, in docking, the minimal energy within the top ten predictions indicates the correct orientation of the construct. Following this criterion helped in designing a better construct.
Constrained peptides such as TDs are good choices for drug development [67]. Smaller cyclic peptides with the same activity of TD have been developed before [68,69], which possessed enhanced antiviral properties and immunosuppressive activity [70].
TDs, as a group of cyclic peptides, have high stability and low toxicity. Moreover, cyclic structures usually show good affinity and selectivity because of the increased surface area and conformational rigidity [68].
In animal models, TDs have demonstrated minimal immunogenicity, low toxicity [35,71], and thermal and proteolytic stability [72]. TDs are β-hairpin antimicrobial peptides, but their amphiphilicity is low, which causes their low hemolytic activity [73]. Compared with alpha-and beta-defensins, TD activities are relatively insensitive to salt, divalent cations, and serum [74]. For instance, a hybrid peptide of TD and human beta-defensin-1 showed saltresistant property [75]. TDs are a type of lectin and can bind to carbohydrate-containing surface molecules, for example, when protecting against HIV-1 infection [76]. Aminoglycosides produce functional TDs in humans [77]. Therefore, they could be used for COVID-19 treatment [78].
One primary reason for COVID-19 severity and mortality is inflammatory cytokine induction at the late stages of the disease, known as cytokine storm or cytokine release syndrome (CRS), which may even lead to organ failure [79]. The organ damage in SARS-CoV-2-infected people seems immune-mediated, similar to autoimmune diseases [80]. Thus, identifying the stage of disease progression and timing of interventions is essential. At the early stages of the disease, it is best to reduce the number of viruses, while hyperinflammation should be controlled at more advanced stages [81,82]. Thus, limiting the inflammatory responses could be as essential as preventing virus replication [83,84] and have a life-saving role in COVID-19 patients [82]. It seems that TD can be helpful in both stages. TD has shown anti-inflammatory activity in vitro and in vivo by inhibiting proinflammatory cytokines' release. In this regard, TD activated the phosphatidylinositol 3-kinase (PI3K)/protein kinase B (AKT) pathway and suppressed proinflammatory signals in immune-stimulated cells. RTD-1 immunomodulation activity was identified in the TLR (Toll-like receptors) and TNF (tumor necrosis factor) pathways via suppressing NF-κB (nuclear factor kappa-lightchain-enhancer of activated B cells) activation and expression as well as cytokine release [85][86][87]. In SARS-infected mice treated with RTD1, decreased levels of proinflammatory cytokines and increased survival were observed [35,36].
Presently, two clinical trials are recruiting patients to study the efficacy of defensins in COVID-19, including (1) a phase II clinical trial on a TD ([97] and (2) a phase II trial on brilacidin [98], a small-molecule human defensin mimetic. The primary mechanism of this molecule is disrupting the viral integrity and entry. Brilacidin has shown a synergistic effect with remdesivir in vitro [99].
Furin inhibitors provoke a lower probability of resistance induction because they do not target the virus [100]. However, furin inhibition might cause toxicity and side effects [101]. Some reports of triggering thrombosis in a mouse model, as the possible adverse effects of defensins, especially alpha-defensin, were reported. Additionally, some defensins promoted viral and bacterial infections in certain biological conditions. Hence, it is vital to consider their probable side effects [101,102]. To reduce their adverse events, pulmonary delivery systems using nanoformulation to deliver antimicrobial peptides [103] could be an option.

Conclusion
Antiviral drugs for COVID-19 are in demand to save the lives of infected people, especially considering various mutations of the SARS-CoV-2 genome. Targeting the host-cell enzymes such as furin can cause less viral resistance. Despite some possible side effects, the advantages of short-term furin inhibition could outweigh [16].
According to our in silico results, TDs can act as furin inhibitors presented by the low electrostatic potential of furin active cleft and high electrostatic potential of the TDs, as well as the appropriate size of TDs for insertion into the furin cleft. Arginine residues have an essential role in TD-furin interaction. This relates to the notion that defensins, as arginine-rich molecules, can inhibit furin. Therefore, TDs could be a potential treatment for COVID-19 by furin inhibition and anti-inflammatory mechanisms. Furthermore, TDs have appropriate properties naturally, such as high stability and low toxicity-however, natural TDs have either high antigenicity or low stability in our analyses. Thus, a novel peptide with optimal antigenicity, stability, electrostatic potential, and binding strength was designed. However, further studies, including in vitro and in vivo experiments, are needed to confirm these results. Additionally, the proposed method could be employed for developing other novel proteins and peptides.

Data Availability
All data generated or analyzed during this study are included in this published article and its Supplementary Information files.

Conflicts of Interest
The authors declare no conflict of interest.

Authors' Contributions
MN was responsible for methodology, investigation, writing, original and final draft preparation, review, and final editing. MRR was responsible for software, data curation, performing the analyses, writing, draft preparation, resources, draft writing, and review and editing. ZM was responsible for the investigation, validation of data, draft writing, formal analysis, writing, and editing. AG was responsible for conceptualization, data curation, visualization, resources, writing, review and editing, supervision, project administration, and funding acquisition. All authors read and approved the final manuscript. Manica Negahdaripour, Mohammad Reza Rahbar, and Zahra Mosalanejad contributed equally to this work.

Acknowledgments
This study has been supported by the Vice-Chancellor of Shiraz University of Medical Sciences under Grant no. 22290.

Supplementary Materials
Supplementary Data 1. Figure S1: cartoon representation of furin cavities. Figure S2: the convexity index for each residue of theta defensins. Figure S3: defensin positive electrostatic potential. Figure S4: a centrality analysis of the important residues involved in the peptide-furin complex. Supplementary Data 2. Table S1: the available UniProtKB data on the furin reference sequence. Table S2: the interactive residues in furin-peptide complexes (available in PepDB). Table S3: the convexity index for each residue of theta defensins. The most protruded side chain is arginine. Table S4: furin-peptide docking, the structure-based method, resulted in ten different orientations for each peptide structure and furin structure. Table S5: designing a novel peptide, which is referred to as "construct." The criteria for designing the construct are achieving a peptide 11 Computational and Mathematical Methods in Medicine with higher stability, lower antigenicity, and higher electrostatic potential. (Supplementary Materials)